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We have implemented the single-site density matrix renormalization group algorithm for the 
variational optimization of SU(2) x U(l) (spin and particle number) invariant matrix product states 
for general spin and particle number symmetric fermionic Hamiltonians. This class also includes non- 
relativistic quantum chemical systems within the Born-Oppenheimer approximation. High-accuracy 
ab-initio finite field results of the longitudinal static polarizabilities and second hyperpolarizabilities 
of one-dimensional hydrogen chains are obtained with the algorithm. A comparison with other 
methods is made. 



I. INTRODUCTION 

The density matrix renormalization group (DMRG) 
outperforms other high-accuracy methods for one- 
dimensional systems. The reason is its underlying ma- 
trix product state (MPS) Ansatz, which efSciently cap- 
tures the quantum entanglement properties of low-lying 
eigenstates. This allows to reach exact diagonalization 
accuracy within a reasonable amount of time. Sections 
In] and inn describe the MPS Ansatz and the DMRG al- 
gorithm. Section IV discusses how global SU(2) x U(l) 
symmetry is imposed on an MPS and how the SU(2) x 
U(l) symmetry of the fermionic Hamiltonian allows to 
work with reduced tensors only. The finite field method 
is explained in sectionjV] The results are shown in section 
|VI| and analyzed in section |VII| Finally, a conclusion is 
made in section |VIII Note that this paper corresponds 
to a poster presentation at the 16th International Con- 
ference on Recent Progress in Many-Body Theories and 
is therefore by no means complete in terms of references. 



II. MATRIX PRODUCT STATES 

A quantum mechanical many-body system consists 
of a Hamiltonian acting on a Hilbert space. The 
Hilbert space is formed as the direct product of L single 
particle spaces. Its dimension increases exponentially 
with L. This has induced the development of many 
approximative solution methods. DMRG is one of them. 
It can be understood in two ways: as a renormalization 
group technique or by means of its underlying Ansatz, 
the MPS. This text focuses on the latter. 

IiP is shown that every quantum state can be rewritten 
as an MPS: 

{Qj} {ij} 

The indices aj represent the local (single particle) de- 



grees of freedom. The indices ij are called the virtual 
indices. For I^E") to represent the full Hilbert space, the 
dimension of the virtual indices has to grow exponen- 
tially toward the middle of the MPS chain. For practical 
calculations, these virtual dimensions are truncated to D 
and the MPS forms a variational Ansatz for the given 
problem. For one-dimensional systems, this Ansatz effi- 
ciently captures the quantum entanglement properties of 
low-lying eigenstates. 



III. THE DENSITY MATRIX 
RENORMALIZATION GROUP 

A general ab-initio quantum chemistry Hamiltonian is 
written in second-quantized form as: 

H ^ Y^{a\f\ PYai^ap, 

+ I {(^P\V\lS)ai^al^asra^, (2) 



where a, j3, 7 and 5 denote orbitals and a and r denote 
spin projections. The global spin and the global particle 
number are conserved by this Hamiltonian. 

The single-site DMRG algorithm sweeps through the 
MPS chain and optimizes at each step a single A-tensor 
variationally, while the others are kept constant. When 
no local minima are encountered, the global minimum of 
(^f I if I VP) is found and the whole MPS is variationally 
optimized. The size of the virtual dimension controls the 
accuracy of the resultP^ 
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Here, a and k are fitting parameters, inexact is the exact 
diagonalization result and Ejj the energy when an Ansatz 
with virtual dimension D is used. The relation ([s]) is 
illustrated in Fig. [T] 
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FIG. 1: The ground state of a hydrogen chain of 36 atoms, 
with an alternate atom spacing of 2/3 a.u., in the Lowdin 
transformed STO-6G basis, is approximated by several MPS 
Ansatze. The scaling of the ground state energy wit h D , 
the virtual dimension per symmetry sector (see section IV I, 
follows 



IV. IMPOSING GLOBAL SU(2) x U(l) 
SYMMETRY 



The use of global symmetries of the Hamiltonian has 
two main advantages^ 

• It enables to scan only the region of interest of the 
many-body Hilbert space. 

• It causes an improvement in algorithmic perfor- 
mance: CPU and memory demands decrease. 

Global symmetry can be imposed by requiring that the 
three-index tensors in the MPS chain are irreducible 



tensor operators of the imposed symmetry grouppBl. The 
local and virtual bases are represented in states with the 
correct symmetry, i.e. spin s or j, spin projection or 
and particle number N. Due to the Wigner-Eckart 
theorem, each A-tcnsor decomposes in a structural part 
and a degeneracy part T: 



(4) 



j3 counts the number of times the irreducible rep- 
resentation {jN) occurs at a virtual boundary^. If 
D{jN) — size(/3), this corresponds to a virtual dimen- 
sion of (2j -I- l)D{jN) in a non-symmetry adapted MPS. 
These virtual dimensions are truncated per symmetry 
sector to D. The desired global symmetry can be 
imposed by requiring that the left virtual index of the 
leftmost tensor in the MPS chain consists of 1 irre- 
ducible representation corresponding to (jlA^l) = (00), 
while the right virtual index of the rightmost tensor 
consists of 1 irreducible representation corresponding to 
UrNr) — {SN), the desired global spin and particle 
number. 



= dl^ and 6™ = (— l)5~'"a_„ trans- 
form as irreducible tensor operators with spin ^ under 
SU(2). All (partial) evaluations of creation/annihilation 
operators can hence be done by implicitly summing over 
the common multiplets and recoupling the local, virtual 
and operator spins. No spin projections or Clebsch- 
Gordan coefficients are used in the algorithrrP2l_ 



V. FINITE FIELD EXTRAPOLATIONS 

The longitudinal static polarizability a^z and second 
hyperpolarizability jzzzz can be obtained for centrosym- 
metric systems asSl 



/2£;(0) - 2E{F)\ 



^ / F-s-O 

-6£;(0) + 8E{F) - 2E{2F)\ 



(6) 
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where E{F) denotes the energy when an electric field F 
in the z-direction is applied. If F is chosen too small, 
noise due to finite precision arithmetic is enhanced. If F 
is chosen too large, higher order effects come into play. 
The minimal finite difference formulas both have a lowest 
order correction term of 0{F'^). We calculate both azz 
and jzzzz for several field values F and extrapolate these 
to F = with 



qiF) = g(0) + cF^ 



(8) 



where q can be ttzz or jzzzz- Values of g(0) and c are 
obtained by the fit. The procedure is shown in Fig. [2] 



VI. RESULTS 

The finite field extrapolation was carried out for hy- 
drogen chains with 

• increasing length: H2N with 5 < A'' < 20 

• different bond length alternations: 2/2.5 a.u., 2/3 
a.u. and 2/4 a.u. 

• the Lowdin transformed minimal basis set ST0-6G 

• 5 levels of theory: HP, MP2, CCSD, CCSD(T) and 
MPS 

The virtual dimension D of the symmetry sectors of the 
MPS was chosen in order that the MPS results match 
the exact diagonalization results. HF, MP2, CCSD and 
CCSD(T) calculations were performed with the molecu- 
lar electronic structure program Dalton'^. The extrapo- 
lation results are shown in Fig. [3j |4] and [5] 
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FIG. 2: Finite field extrapolations for MPS calculations of a hydrogen chain of 36 atoms, with an alternate atom spacing of 
2/3 a.u., in the Lowdin transformed ST0-6G basis. 
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FIG. 3: Polarizabilities and second hyperpolarizabilities of hydrogen chains with an alternate atom spacing of 2/2.5 a.u., in 
the Lowdin transformed STO-6G basis, for different levels of theory. 



VII. DISCUSSION 

The change in bond length alternation results in a 
different separation of H2 molecules, which themselves 
have a fixed bond length of 2 a.u. The multireference 
character doesn't change significantly for the three bond 
length alternations, but the correlation length of the 
electrons in the lattice does. 

The results for a^z and jzzzz are ordered in magnitude 
as: 

CCSD < CCSD{T) < MPS < MP2 < HF (9) 

This ordering differs from^, where an aug-cc-pVDZ ba- 
sis set was used for the HF, MP2, CCSD and CCSD(T) 
methods. The ordering is hence basis set dependent. The 
results for azz and ^zzzz are ordered according to their 
deviation from the exact diagonalization results as: 

HF > MP2 > CCSD > CCSD{T) > MPS (10) 

which confirms that CCSD(T) is more accurate to cal- 
culate the polarizability and second hyperpolarizability 



than HF, MP2 and CCSE0. 

The deviation between the CCSD(T) and exact diago- 
nalization results for azz and 'jzzzz is the largest for the 
2/2.5 a.u. bond length alternation. The electrons in the 
2/2.5 a.u. configuration are more mobile throughout the 
chain than the electrons in the 2/4 a.u. configuration. 
They are thus correlated on a longer length scale in the 
lattice, i.e. with more electrons. Results obtained with 
a method that captures correlation only partially are ex- 
pected to deviate more from the true results in that case. 
Fortunately, the amount of electron correlation captured 
by the MPS Ansatz is controlled by the truncation of the 
virtual dimension and can be finetuned to any desired 
accuracy. 



VIII. CONCLUSION 

The MPS Ansatz is able to capture all required corre- 
lation to be able to obtain highly accurate ab-initio finite 
field results of the longitudinal static polarizability and 
second hyperpolarizability of one-dimensional chains. 



4 



14 
13 
12 

9 



2/3 a.u. 



f.04 0.06 0.08 0.10 0.12 0.14 0.16 0.18 0.20 0.22 
(N = number of JJa-dimers}"^ 




r04 0.06 0.08 0.10 0.12 0.14 0.16 0.18 0.20 0.22 
{N = number of /f.,-dimers)"^ 



FIG. 4: Polarizabilities and second hyperpolarizabilities of hydrogen chains with an alternate atom spacing of 2/3 a.u., in the 
Lowdin transformed ST0-6G basis, for different levels of theory. The legend is shown in Fig. |3] 
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FIG. 5: Polarizabilities and second hyperpolarizabilities of hydrogen chains with an alternate atom spacing of 2/4 a.u., in the 
Lowdin transformed ST0-6G basis, for different levels of theory. The legend is shown in Fig. |3] 
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itly kept in the reduced basis to avoid losing quantum 
numberJ^. Note that the division of the virtual indices in 
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the Wigner-Eckart theorem to work with reduced tensors. 



